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ABSTRACT 

We estimate the impact of the Dipole Strayhght Contamination (DSC) for the 
Planck satelhte on the ahgnnients of vectors associated to the low multipoles 
of the pattern of the cosmic microwave background (CMB) anisotropies. In 
particular we study how the probability distributions of eighteen estimators for 
the alignments change when DSC is taken into account. We find that possible 
residual DSC should leave a non-negligible impact on low multipole align- 
ments for effective values of the fractional far sidelobe integrated response, p, 
larger than ~ few x 10"'^. The effect is strongly dependent on the intrinsic 
sky amplitude and weakly dependent on the considered scanning strategy. We 
find a decrease of the alignment probability between the quadrupole and the 
dipole and an increase of the alignment probability between the hexadecapole 
and the dipole (larger is the intrinsic sky amplitude and lower is the contam- 
ination). The remaining estimators do not exhibit clear signatures, except, in 
some cases, considering the largest values of p and the lowest sky amplitudes. 
Provided that the real sidelobes of the Planck receivers in flight conditions 
will correspond to p <, few x 10"'^, as realistically expected at least in the cos- 

mological frequency channels, and will be known with accuracies better than 
~ few X 10% allowing for a suitable cleaning during data reduction, Planck 
will be very weakly affected from DSC on the alignments of low multipoles. 

Key words: Cosmology: cosmic microwave background - methods: data 
analysis. 



1 INTRODUCTION 

The anisotropy pattern of the Cosmic Microwave 
Background (CMB), obtained by Wilkinson Microwave 
Anisotropy Probe (WMAF|3), probes cosmological 
mode ls with unprecedented precision IjSpergel et al.l 
I2OO6I '). Although WMAP data are largely consistent 
with the concordance A Cold Dark Matter (ACDM) 
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model, there are some interesting deviations from it, 
in particular on the largest angular scales: a surpris- 
ingly low amplitude of the quadrupole term of the 
angular power spectrum (APS), found by Cosmic 
Background Explor er (COBE) (jSmoot et al' '1990!; 
iHinshaw et"al]|l996l ) and WMAP fSpcrgcl ct al. 2006), 
and an unlikely (for a statistically isotropic random 
field) alignment of the quadrupole and th e oc tupole 
(ITegmark. de Oliveira Costa and Hamiltor] 20031: 



2005a: 



Copi et al.ll2004l : ISchwarz et al."2004ULa nd fc Magueiiol 



Vale I 



2005: 



Abramo et al.l 



20061 ). 



Weeks 2005; Co pi et al.l I2OO6I : 
Moreover, both quadrupole and 
octupole align with the CMB dipole. Other unlikely 
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alignments are present in the aforemen tioned papers 
and ot her low I anomalies are described in lEriksen et alj 
(|2004l ). 

It is still unknown if these anomalies come from fun- 
damental physics or if they are the residual of some not 
removed astrophysical or systematic effect. This open 
question has attracted a lot of interest and many pa- 
pers have been published about this subject in the last 
few years. 

T his paper represents a step forward o f a previous 
work (iBurigana. Gruppuso fc Finellil l|2006l ). henceforth 
BGF06) where the impact of the systematic effect in- 
duced by the CMB kinematic dipole signal entering the 
main spillover (Dipole Straylight Contamination, DSC) 
on the APS has been studied in particular for the forth- 
commg Planc^ mission. If ere we wish to estimate the 
main implications of the same systematic effect on the 
issue of alignments of the low multipole vectors under 
experimental conditions (observational strategy, main 
properties of the far sidelobes) like those typically pre- 
dicted for Planck in its cosmological frequency channels 
(at the lowest and highest Planck frequency channels 
Galactic straylight dominates over dipole straylight). 

The measurement of the low I pattern is affected 
by cosmic variance, foregrounds and systematic effects 
[for a d iscussion on Ga lactic straylight c ontamination 
see e. g. iBurigana et al.l (|200I, 20o"3 ) and ISandri et al.l 
(|2004l ) in the co ntext of the Planck Lo w Frequency In- 
strum ent (LFI) iMandolesi et al.]|l998l ). iLamarre et al.l 
(|2002 ) in the con text of the Planck High Frequency In- 
strum e nt (HFI) (Puget e t al. 1998), and iBarnes et al.l 



for the dipole. 



(Hool) , IChiang. Coles fc Naselsky [2QQ(h, in the context 



of WMAP; see e.j 



IChiaiig. Naselskv fc Coles! (|2006D : 
iNaselskv fc Verkhodanovl (120061 ) for an analysis of large 
angular scale foreground contamination in WMAP 
data]. 

Through a simple analytical model [top-hat 
approximation for the ma i n spi llover response 
ijGruppuso, Burigana fc Finellil 120061 )] and several 
numerical simulations we tackled the systematic effect 
induced on the APS at low and intermediate multipoles 
by the CMB kinematic dipole signal entering the main 
spillover (BGF06). In that study, we analytically found 
that in one survey or in an odd number of surveys 
the DSC map, described by the coefficients af^, that 
turn on for ^ < 4, is given by 



SL 




C\Ol , 



SL 1 StT 



(1.1) 



(1.2) 



^ http://www.rssd.esa.int/planck 

^ With "number of surveys" we mean "number of full sky 
mappings" consecutively realized by the satellite. 



2id2) C23 : 



for the quadrupole and 



SL 



SL 
0-4±4 




± 2id2) C23 , 



12 /35 , , , , , 

\ ("1 ± 4ja2 ) C23 . 

225 V 27r ^ ^ ■ 



for the hexadecapole, where C23 = C2 + C3 and 



ci = ^/y^Tv/sLA sin(2A) Tio , 

C2=4y^3/8^fsLA, 

C3 ^ 4^/3/8^ fsL sin(2A)/2, 
di = sin A Tie [Til] , 
d2 = sin A 7m [Tn] , 



(1.3) 



(1.4) 



(1.5) 



(1.6) 
(1.7) 
(1.8) 
(1.9) 
(1.10) 



with fsL = p/(4Asin A), p being the ratio between the 
power entering the main spillover and the total power 
entering the receiver (i.e. essentially the power entering 
the main beam) and A being the angular side of the 
box that in the {8, <^)-plane describe the main spillover 
region. Note that the octupole is unaffected. The terms 
Tio, Re [Til] and Im [Tii] are the coefflcients of the kine- 
matic dipole. In ecliptic coordinates their values in ther- 
modynamic temperature are 



Im [Til] = 0.69823 mK, 
Tio = -1.32225 mK, 
Re[Tii] = 4.69963 mK. 



(1.11) 
(1.12) 
(1.13) 



The expressions for af^ are obtained perturbatively to 
the first order in the angle (a) between the directions of 
the main spillover and of the spin axis. Moreover, it has 
been supposed that the main spillover centre is located 
on the plane defined by the spin axis and the telescope 
line of sight {0mb,'Pmb)- The relaxation of the latter as- 
sumption is described by the introduction of a phase /3 
that parametrize the displacement of the main spillover 
direction from 9mb- The above expressions hold in the 
case of a simple scanning strategy with the spacecraft 
spin axis always on the ecliptic plane, i.e., in the case 
of Pla nck, for the so-called no minal scanning strategy 
(NSS) (|Dupac fc Taubeij|2005l '). The treatment of com- 
plex scanning strategies require numerical simulations. 

Note the simple pattern at low £ due to DSC: 
even multipoles are modified at the leading order in 
a whereas odd multipoles do not change or are only 
weakly contaminated (i.e. linearly in a). The introduc- 
tion of /3 does not change this scheme and /3 appears 
only to the linear order in a in the odd multipoles 
(BGF06). Therefore, for the dipole the introduction of 
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/3 leads to the replacements 



SL SL o 1 ^ / 3 J ■ r, 

flio ^0.10 cos/i + -W — d2C23asin^, 
3 V 47r 



SL SL o , 16 / 3 ■ ■ r, 

3 V Stt 



(1.14) 



(1.15) 



In addition, for the octupole we have the following non 
vanishing coefficients 

2 / 14 

af±2 = TT- (±2idi + d2) C23Qsin/3 , (1-16) 

3 V tOTT 



St 16 / 7 . . 
a3±3 = TTA/ ^*ciQsm/3. 

9 V OTT 



(1.17) 



In order to investigate the implications of a non- 
proper subtraction of DSC on the low I alignments for 
the Planck experiment, we consider not only the above 
summarized analytical description (where the NSS is 
adopted) but also numerical simulations to explore the 
case of a cycloidal scanning strategy (CSS) (with slow 
precessions) that is beyond the analytical approxima- 
tion. 

We shall consider the case of a single survey (or of 
an odd number of surveys), so providing upper limits 
to the contamination induced by this effect. Except for 
the dipole, the final DSC impact is in fact significantly 
reduced by considering an even number of complete sur- 
veys, although in a way dependent on the considered 
scanning strategy (we remember that also in the case of 
an even number of surveys a remarkable effect survives 
to the averaging if one survey is not complete). 

This article is organized as follows. In Section [5] the 
multipole vectors expansion and its link with spherical 
ha rmonic e xpans ion is reviewed along the line suggested 
by IWeeksl (|2005l 1^ discussing also the complementarity 
between the information contained in the APS and in 
the multipole vectors. In Section[3]we describe the simu- 
lations that we have performed and how the estimators 
for the statistical analysis are defined. In Section |4] we 
present how the probability distribution of the estima- 
tors change when DSC is taken into account. Finally, we 
discuss the obtained results and draw our conclusions in 
Section (5] 



2 MULTIPOLE VECTORS 

The alignment of multipoles is better understood by a 
new representation of CMB anisotropy maps where the 
atm (coefficients of the expansion over the ba sis of spher- 
ical harmonics) are replaced by vectors (|Copi et al.l 



http: / /www. geometrygames.org/Maxwell/ 



120041 '). In particular, each multipole order I is repre- 
sented by I unit vectors and one amplitude A 



A 



Ml, 



(2.18) 



Note that the number of independent objects is the 
same in the l.h.s and r.h.s. of equation (|2.18p : 2£ + 1 
for atm equals 31 (numbers of components of the vec- 
tors) 4-1 (given by A^^^) —£ (because there are £ con- 
straints due to the normalization conditions of the vec- 
tors). One of the advantage of this representation is that 
from these unit vectors one can easily construct scalar 
quantities that are invariant under rotation. Note that 
is not equally easy to obtain scalar quantities directly 
from the aim coefficients that, of course, depend on the 
coordinate system. 

Equation (|2.18ll can be understood starting from 
this observation (Wc oka. 2005. 1: if / is a solution of the 
Laplace equation 



(2.19) 



dx +dy+d2 in Cartesian coordinates, then 



where ^d^+dy 
it is possible to build a new solution /' applying a di 
rectional derivative to / 

V^/' = , (2.20) 



Va/ = u- V/ = /', 



with the gradient V = {dx,dy,dz). This hap pens be- 
cause the two operators and Va commute. iMaxwelll 
(|l873l ') repeated this observation £ times considering the 
1/r potential as starting solution. Here r = {x,y,z) and 
r = Vr-f — ^ -\- ip- -\- . In this way, one obtains 



fe{x,y,z) = Vu/, ■ ■ ■ Va^Vsi - . 



(2.21) 



Observe the simple pattern that emerges as we apply 
the directional derivatives one at a time: 



,fo 
h 
h 
h 



(-l)(^i-^") 

(3 ■ l)(jii ■ r){u2 ■ r) + r^(— Ul ■ U2) 
(-5 ■ 3 • l)(ui • r){u2 ■ r){u3 ■ r) + r-2(...) 



The (...) stands for a polynomial which we do not write 
explicitly, being useless for the current purposes. 

Moreover, writing ft in spherical coordinates once 
r is set to 1, one finds the following property 



\/^Mi,e,(b) = £i£ + i)Mi, 



(2.22) 



where V'^ is the angular Laplace operator defined as 



^i- de (smOde) 
smO 



+ ■ 



■di 



(2.23) 



In other words /^(l, 6, (j)) is eigenfunction of the angular 
part of the Laplace operator with eigenvalue given by 
£(£ 1). This is nothing but the definition of spherical 
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harmonics Yi „ (|Sakurai|[l985l) . Therefore, for every f 
we can write 

= a,™F^„(e,0), (2.24) 

m = — £ 

where the amphtude yl'*' has been inserted because of 
normahzation. Equation (|2.24p makes evident the as- 
sociation represented by equation (|2.18|l . From equa- 
tion (|2.24[) it is possible to write down the set of equa- 
tions that has to be solved to pass from airn to mul- 
tipole vectors. In order to see that this set is solvable 
we count the equations and the unknowns involved in 
this set. From equation (|2.24p we have 21 + 1 equations 
(one equation for each independent aim El) plus t equa- 
tions from the normality conditions of the vectors (i.e. 
Ui-Ui — 1 where i runs from 1 to ^). Therefore the total 
number of independent equations is 3^-1-1. This is also 
the number of unknowns because we have 3 unknowns 
for each vector plus 1 given by the amplitude A'^' . This 
shows that the set is solvable. 

Unfortunately, an explicit analytical solution is pos- 
sible only for I — 1 and for £ ^ 1 numerical method are 
needed 0. For ^ = 1 we have 
1 

-A'''\d-r) = J2 ciimY^^{e,<l>), (2.25) 

m— — 1 

d-d = 1. (2.26) 

Considering the expression for Yim{e, 0) (|Sakurailll985l ) 
we have the following analytical solution 



da: = 




, (2.27) 


dy = 


±«iiV\/«?o/2 + ((4f)^ + («ff)^), 


(2.28) 


dz — 




(2.29) 






(2.30) 


where d = 


{dx,dy,dz) and the labels (R) and (/) stand 



for real and imaginary part. 

An elegant way of solving numerically the above 
set of equations (in order to obtain the multipole vec- 
tors expansion from the set of aem ) is p resented in 
iKatz fc Week3 (|2004| ) and in IWeekd Hooi) where the 
problem of finding £ vectors is translated into the prob- 
lem of finding the zeros of a polynomial of degree 2£. 

^ In fact we would have 4£ -|- 1 equation because each £ dif- 
ferent from has a real and imaginary part. But consider- 
ing that aim with m > are related to those with m < 0) 
through a*^ = {—l)"^ai_m we are left with 2£+l equations. 
® Indeed, for £ = 2 it is possible to obtain the multipoles vec- 
tors computing the eigenvectors of a sym metric and traceless 
tensor represen t ing t he quadrupole [see iLand &: Magueiiol 
l l2005bl ') lDennii ||2005| ')]. 



T his meth od has been implemented in a code developed 
bv lWeeksl (HoOS) whose use is aknowledged. 
We end this section with two observations: 

• equation (|2.24|) is invariant under the change of 
sign of an even number of Ui or under the change of 
sign of an odd number o f Ui and A^^^ (jAbramo et al.l 
l2006l : lKatz fc WeeksHiooi ). This "reflection symmetry" 
implies that in fact multipole vectors define only di- 
rections. Of course the same symmetry has to be satis- 
fied by quantities defined through the multipole vectors. 
Hence the estimators introduced in Section |3l are sensi- 
tive only to directions of multipole vectors (or in other 
words, they are invariant under change of the sign of 
the vectors) . Notice the ambiguity of sign for the dipole 
in equations (|2.27l2.30p as an example of the above 
mentioned refiection symmetry. 

• equation (|2.24l) is also invariant under the trasfor- 
mation 

aim caim , (2-31) 
^W^^c^W^ (2.32) 

where c is a constant. Therefore the £ multipole vec- 
tors associated to aim are the same as those that are 
associated to caim- This means that for a sufficiently 
large number N of random extractions of the 2£ + 1 in- 
dependent values of aim, the corresponding N sets of £ 
multipole vectors do not depend on Ci, the variance of 
aim, in the remarkable case in which the aim follow a 
Gaussian distributional. As a consequence, the same ap- 
plies to the estimators for the alignments (see Section|3}. 
For instance if aim follow a Gaussian distribution, the 
information contained in the multipole vectors is "or- 
thogonal" to that contained in the APS. This property 
is broken in the presence of a systematic effect alter- 
ing the above symmetry. In this case the global effect 
could depend on the intrinsic value of Cf'^ . In partic- 
ular, if not fully removed, DSC will provide a spurious 
deviation from Gaussianity. 



3 ANALYSIS 

In order to study the effect of DSC on the alignment 
of low multipole vectors we have extracted 3 x 10^ 
sky realizations for two different APS amplitudes cor- 
responding to a concordance ACDM model and to 
WMAP. In other words, we have extracted for 
^ = 2, ^ = 3 and I = A from a Gaussian distribu- 
tion with zero mean and variance given by C2, C3, 
C4 and each extraction has been contaminated by af^. 
In particular defining ATi = £{£ -I- l)Ce/{2-K) we have 
taken AT2 = 1250 ^K^ AT3 = 1150 ^tK^ and AT4 = 

This is true for every distribution writable as j ifliml 
where cr is a parameter. 
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Figure 1. Uncontaminated distribution for all the considered estimators. In order from left to right and above to below: W2, 
W2, W2°, R22, W3, W3, D23, S23, D33, 533, Wi, Wi, D42, 542, Z)43, 543, 1)44, 544. All the panels present the counts 
(j/-axis) versus the statistic (x-axis). Vertical lines stand for WMAP value reported in Table 1. For R22 we overplot the analytic 
distribution. See also the text. 



1110 ^iK^ for the ACDM case and AT2 = 211^K^ 
ATs = 1041 /iK^ and AT4, = 731 /xK^ for the WMAP- 
Uke ampUtude case l|Hinshaw et all I2OO6I '). Moreover, 
we set /3 = for sake of simphcity and considered 
p = 1/1000, 4/1000, 7/1000, 10/1000, i.e. in a range rep- 
resentative of current CMB anisotropy space experi- 
ments Uke WMAP and Planck (see e.g. | Barnes ct al 
(|2003l ') and ISandri et all Hooi)). The other parame- 
ters are freezed to the following values a = 7r/18, 

© 0000 RAS, MNRAS 000, 000-000 



A — Tv/W, for sake of simplicity. The coefficients a'^^ 
and 4- af^ are then transformed to multipole vec- 
tors through the Weeks' code. 

As written above, not only the analytical expres- 
sions for af^ have been used to describe DSC but also 
numerical simulations have been performed in order to 
consider cases beyond the analytical approximations. In 
particular, we have numerically taken into account the 
case of a CSS with slow precessions as described for 
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Estimator 


WMAP value 


Probability % 


Estimator 


WMAP value 


Probability % 


W2 


0.901 


9.96 


533 


0.480 


21.76 


W2° 


25.68° 


9.94 


W4 


0.528 


38.73 


W2 


0.854 


8.33 


W4 


0.479 


36.01 


R22 


0.319 


48.04 


DA2 


0.422 


76.82 


W3 


0.796 


3.92 


542 


0.373 


62.80 


W3 


0.735 


2.03 


D43 


0.447 


86.43 


D23 


0.829 


2.74 


543 


0.384 


69.75 


523 


0.732 


0.96 


D44 


0.397 


82.03 


D33 


0.554 


27.49 


544 


0.330 


68.70 



Table 1. WMAP values for the considered estimators adopting the CMB component in the V band of the WMAP 3 year 
release (see the text for further detailes). The percentages are the probability that a random map would have values for the 
estimators larger than those observed in the V band by WMAP. Only for 14^2° the percentage represents the probability to 
extract randomly a value lower than that obtained from the above mentioned data set. 



example in lDupac fc Tauberl (|2005l ). assuming a period 
T = 6 months and a semi-amplitude of 10° , and consid- 
ering a — 10° and /3 = 60°. In this case, we parametrize 
the main spillover response according to the Gaussian 
approximation. With reference to § 3.3 of BGF06, we 
note that in the remarkable case of pencil beam the top- 
hat and Gaussian approximations are essentially equiva- 
lent for small A. Differently, the two cases are not equiv- 
alent for general values of a. On the other hand, it easily 
to show that in the case a = they are equivalent up to 
second order in a and A provided that = A^/3 (and 
b = p). We then adopt here this relation to define the 
beamwidth, a, of the main spillover in the Gaussian ap- 
proximation given the chosen value of A in the top-hat 
approximation used in the analytical approach. 

We present in Section|3]how the distributions of the 
eighteen estimators change when the DSC is properly 
taken into account. The considered estimators are: for 
the alignment quadrupole-dipole 



W2 = \q-d\, 
W2 = \q-d\, 
W2° = arccos(|g-d|)180/7r, 

for the self alignment of the quadrupole 

i?22 = |52i -9221, 

for the alignment octupole-dipole 

3 

W3 = ^|o>-d|/3, 

i = l 
3 



W3 = ^|o, -dl/S, 



(3.1) 
(3.2) 
(3.3) 

(3.4) 

(3.5) 
(3.6) 



for the alignment quadrupole-octupole 

3 

S23 = ^|g-o,|/3. 



i=l 
3 



D23 = ^!<7- 6,1/3, 
1=1 

for the self-alignment of the octupole 

3 

5*33 = \oi-Oj\/3, 

3 

-D33 = \6,-dj\/3, 

for the alignment hexadecapole-dipole 

6 

i = l 
6 

WA = ^je, -dj/e, 

i = l 

for the alignment hexadecapole-quadrupole 

6 

^42 = ^|e,-g-l/6, 
1=1 

6 



D42 = ^|e,-qj/6, 



for the alignment hexadecapole-octupole 

6 3 

i=i j=i 



(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 
(3.12) 

(3.13) 
(3.14) 

(3.15) 
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6 3 

= ^ 51 -^^1/18' (3.16) 

i=i j=i 

and for the self- alignment of the hexadecapole 

6 

S44 = ^ IS-e-jl/lS, (3.17) 

i — 1 ,j >z 
6 

_D44 = ^ \e,-ej\/15, (3.18) 

i — 1 ,j >z 

where the symbol "hat" stands for a vector with norm 
equal to 1 and where the "area vectors" are defined as 



hi X 922 , 


(3.19) 


Ol — 632 X 633 , 


(3.20) 


02 = 633 X 631 , 


(3.21) 


03 = 631 X 632 , 


(3.22) 


ei = (=41 X 642 , 


(3.23) 


6*2 = 641 X 643 , 


(3.24) 


6*3 = 641 X 644 , 


(3.25) 


6*4 = 642 X 643 , 


(3.26) 


6*5 = 642 X 644 , 


(3.27) 


6*6 — 643 X 644 , 


(3.28) 



with q2j representing the two normahzed multipole vec- 
tors (j = 1, 2) associated to the quadrupole, 63^ rep- 
resenting the three normalized multipole vectors (j = 
1,2,3) associated to the octupole and 64-, representing 
the four normahzed multipole vectors (j = 1, 2, 3, 4) as- 
sociated to the hexadecapole. Notice that all the estima- 
tors but W2° G [0°, 90°], belong to the interval [0, 1] and 
contain absolute values in order to make them invariant 
under the refiection symmetry discussed in Section [21 
Notice also that i?22 is the unique es timator whose dis- 
tribution is analytica lly known [see iLand fc Magueiiol 
(|2005bl ) lDenni3 (200511 ]: 

where x = i?22. 



4 RESULTS 

In Fig. [1] we plot the uncontaminated distribution for 
W2, W2, W2°, R22, W3, W3, D23, S23, D33, 533, 
#4, W-i, D42, 542, D43, 543, D44, 544. Notice the 
perfect agreement between the analytical, see equa- 
tion (|3.29p . and the numerical distribution of _R22. For 
sake of completeness we report in Table [T] the values of 
the considered estimators obtained, for example, adopt- 
ing the CMB anisotropy component in the V band of 
the WMAP 3 year rele ase wi t h a K p2 mask and the 
dipole given in iHinshaw et al.l l|2006l ). 



In all the subsequent figures we display the differ- 
ence (DD) between the distribution contaminated by 
DSC and the uncontaminated one: Figs 1 ^141 refer to 
the considered NSS. The Figures for the considered CSS 
are not shown being very similar to the NSS. However 
some (small) differences are reported in Section [S] Of 
course, the sum of counts of the uncontaminated distri- 
bution is always 3 x 10* while the sum of counts of the 
DD is zero. 

The DD for W2, W2, W2° are shown in Fig. [2] for 
the considered ACDM model and in Fig.|3]for the model 
with WMAP amplitude. For W'i, W3, L>33 and 533 the 
DD is exactly zero for the considered NSS. Therefore 
they are not reported for sake of conciseness. For the 
considered CSS they are noisy-like. The DD for _R22 are 
plotted in Fig|4]for the considered ACDM model and for 
the WMAP amplitude. The DD for D23 and 523 are 
shown in Fig. [S] for the considered ACDM model and 
in Fig. [6] for the WMAP amplitude. The DD for W4: 
and WA: are plotted in Fig. [7] for the considered ACDM 
model and in Fig. [S] for the WMAP amplitude. The DD 
for _D42 and 542 are given in Fig. [9] for the considered 
ACDM model and in Fig. [10] for the WMAP amplitude. 
The DD for D43 and 543 are shown in Fig. [TT] for the 
considered ACDM model and in Fig. [12] for th WMAP 
amplitude. The DD for _D44 and 544 are plotted in Fig. 
[T^ for the considered ACDM model and in Fig. [2] for 
the WMAP amplitude. The not reported DD for NSS 
are essentially noisy-like. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have studied the impact of a non- proper 
subtraction of DSC on alignments of multipole vectors 
associated to low multipoles for the forthcoming Planck 
mission. This work represents a step forward of BGF06 
in the study of DSC on the low multipoles of CMB pat- 
tern. 

The analyses presented in this paper and in BGF06 
provide a working example of a possible connection be- 
tween the two low I anomalies (APS and alignments 
of the low multipole vectors) in the presence of a non 
cosmological residual in the data. To our knowledge, 
a non-properly subtracted systematic effect (of instru- 
mental or astrophysical origin, or, as in this case, com- 
ing from a combination of them) represents the easiest 
way to link the statistics of aim and the amplitude of 
Ci, otherwise disconnected since the symmetry defined 
in equations (|2.31|l and (|2.32|l holding at least for any 
distribution writable as in footnote 7 (like the Gaussian 
one). 

We summarize the main results of this paper sepa- 
rately for the considered nominal and cycloidal scan- 
ning strategies in the context of the Planck mission 
(|Dupac fc TaubeJl2005l '). 

For the NSS: 
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Figure 2. Differences between the DSC contaminated distribution and the uncontaminated one in ACDM model for the 
alignment Quadrupole-Dipole in the case of small a, vanishing f) and nominal scanning strategy. First row: W2. Second row: 
W2. Third row: W2° . From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels 
present the counts (j;-axis) versus the statistic (x-axis). See also the text. 




Figure 3. Differences between the DSC contaminated distribution and the uncontaminated one for the WMAP amplitude for 
the alignment Quadrupole-Dipole in the case of small o, vanishing (3 and nominal scanning strategy. First row: 14^2. Second 
row: W2. Third row: W2° . From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels 
present the counts (y-axis) versus the statistic (x-axis). See also the text. 



• the probability of alignment quadrupole-dipole 
(#2,1^2,^2°) tends to be lowered; the impact of this 
effect increases going from ACDM-like APS amplitudes 
to WMAP-like APS amplitudes and for increasing val- 
ues of p, the percentage of power entering the main 
spillover with respect to the main beam (see Figs [2] and 

ED; 

• some features show up in the estimator for the 
alignment hexadecapole-dipole (M^4,W^4) if p is suffi- 



ciently large, both for ACDM and WMAP-like ampli- 
tudes (see Figs [7] and [S| ; 

• some features show up in the estimator for 
the alignment hexadecapole-quadrupole (D42,S'42) but 
only for large p, i.e. p = 1/100, and WMAP-like intrinsic 
amplitude (see Fig. I10|l ; 

• the remaining estimators {R22, D23, S23, D42, 
542, D43, 543, D44, 544) do not show remarkable fea- 
tures, being essentially noisy-like (see Figs |5] [9] 1111 1121 
1131 and I14p : on the other hand, a weak signature ap- 
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Figure 4. Differences between the DSC contaminated distribution and the uncontaminated one for the self-alignment 
Quadrupole-Quadrupole in the case of small a, vanishing /3 and nominal scanning strategy. First row: R22 in ACDM model. 
Second row: R22 for the WMAP amplitude of the intrinsic sky. Prom left to right (in every row) p = 1/1000, p = 4/1000, 
p = 7/1000, p = 10/1000. All the panels present the counts (j/-axis) versus the statistic (x-axis). See also the text. 




Figure 5. Differences between the DSC contaminated distribution and the uncontaminated one in ACDM model for the 
alignment Octupolc-Quadrupole in the case of small a, vanishing /3 and nominal scanning strategy. First row: D23. Second row: 
523. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present the counts 
(y-axis) versus the statistic (x-axis). See also the text. 



pears in the case of R22, S23 and S'44 (for the alignment 
octupole-quadrupole and self- alignment of the hexade- 
capole, respectively) for WMAP-like intrinsic ampli- 
tudes and the maximum considered value of p (see Figs 
Il[6]and[l4|. 

For the CSS: 

• the alignment quadrupole-dipole {W2,W2,W2°) 
again tends to be lowered, the impact of this effect be- 
ing again stronger if the intrinsic sky amplitude is lower 
and p is sufficiently large ; 

• some features appear for the estimator of the align- 
ment hexadecapole-dipole (1^4,^4) but stronger than 
those obtained for the NSS and increasing with p ; 

• there is a clear signature in the estimator of 
the alignment hexadecapole-quadrupole (D42,5'42) but 
only for WMAP-like amplitudes ; 

• again, the remaining estimators {R22, D23, S23, 
D43, 543, D44, 544) do not show remarkable fea- 
tures, being essentially noisy-like; on the other hand, 
a weak signature appears also for this scanning strat- 
egy in the case of R22, S23 and S44 (self alignment 



of the quadrupole, alignment octupole-quadrupole and 
self alignment of the hexadecapole respectively), for 
WMAP-like intrinsic amplitudes and the maximum con- 
sidered value of p. 

We conclude that possible residual DSC should 
leave a non- negligible impact on low multipole align- 
ments for far sidelobe integrated responses correspond- 
ing to effective values of p ^ few x 10"'^. 

We note also that, in general, it could be very 
useful to carry out the alignment analysis by exploit- 
ing both normalized and unnormalized estimators, since 
they could show different sensitivity in the diagnostic of 
the same effect, as evident for example in Fig. (5] 

The scanning strategy and the behaviour of the op- 
tics of WMAP are significantly different from those con- 
sidered here, suitable for Planck. It is then difficult to 
extrapolate the above results to the WMAP surveys. 
If the large scale footprint of the DSC is not so criti- 
cally dependent on the above details and our analysis 
applies as well to the WMAP surveys, DSC may not 
easily explain the whole set of anomalous alignments 
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Figure 6. Differences between the DSC contaminated distribution and the uncontaminated one for the WMAP amplitude for 
the ahgnment Octupole-Quadrupole in the case of small a, vanishing f3 and nominal scanning strategy. First row: D23. Second 
row; S23. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present the counts 
(y-axis) versus the statistic (x-axis). See also the text. 
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Figure 7. Differences between the DSC contaminated distribution and the uncontaminated one in ACDM model for the 
alignment Hexadecapole-Dipole in the case of small a, vanishing /3 and nominal scanning strategy. First row: W4. Second row: 
W4. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present the counts 
(j/-axis) versus the statistic (x-axis). See also the text. 



at l arse scale found a lso in WMAP three year release 
(see ICopi et al] (|2006l )). However, we believe that fur- 
ther studies of this and other systematic effects are re- 
quired. 

Provided that the real sidelobes of the Planck re- 
ceivers in flight conditions will correspon d to values of 
p ^ few X 10~^ - as realistically expected l|Sandri et al.l 

I2OO4I ') at least in the cosmological frequency channels - 
and will be known with relative accuracies better than 
~ few X 10% (leaving to smaller residual contamina- 
tions, equivalent to p ^ 10"'^, after a suitable cleaning 

during data reduction) Planck maps will be very weakly 
affected by DSC on the alignments. 
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Figure 10. Differences between tlie DSC contaminated distribution and tlie uncontaminated one for the WMAP amplitude 
for the alignment Hexadecapole-Quadrupole in the case of small a, vanishing /3 and nominal scanning strategy. First row: £'42. 
Second row: S'42. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present 
the counts (j/-axis) versus the statistic (x-axis). See also the text. 
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Figure 11. Differences between the DSC contaminated distribution and the uncontaminated one in ACDM model for the 
alignment Hexadecapole-Octupole in the case of small a, vanishing (3 and nominal scanning strategy. First row: DAS. Second 
row: 543. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present the counts 
(y-axis) versus the statistic (x-axis). See also the text. 
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Figure 12. Differences between the DSC contaminated distribution and the uncontaminated one for the WMAP amplitude 
for the aUgnment Hexadecapole-Octupole in the case of small a, vanishing /3 and nominal scanning strategy. First row: D43. 
Second row: 543. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present 
the counts (j/-axis) versus the statistic (a;-axis). See also the text. 
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Figure 13. Differences between the DSC contaminated distribution and the uncontaminated one in ACDM model for the 
self-alignment of the Hexadecapole in the case of small a, vanishing /3 and nominal scanning strategy. First row: Z)44. Second 
row: S44. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present the counts 
(y-axis) versus the statistic (x-axis). See also the text. 
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Figure 14. Differences between the DSC contaminated distribution and the uncontaminated one for the WMAP amplitude 
for the self-alignment of the Hexadecapole in the case of small a, vanishing /3 and nominal scanning strategy. First row: D44. 
Second row: S44. From left to right (in every row) p = 1/1000, p = 4/1000, p = 7/1000, p = 10/1000. All the panels present 
the counts (y-axis) versus the statistic (a;-axis). See also the text. 
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